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We used neutron powder diffraction to obtain the local and long-range structure of 
Lao.sCao.sMnOg at room temperature and 20 K. By combining Rietveld and pair distribution func- 
tion analysis of the total neutron scattering data, we have analyzed the structure of the compound 
using two competing models describing the low temperature phase: first the charge ordered/orbital 
ordered model and second the Mn-Mn dimer model. These structural models fit the 20 K neu- 
tron powder diffraction pattern equally well using Rietveld analysis. Therefore, pair distribution 
function analysis is used to probe the local and medium-range structure revealing a system with 
two distinctly distorted Mn octahedra and Mn ions with non-integral valence states. The distorted 
octahedra differ with the structural model for the Zener polaron type Mn-Mn dimer picture pro- 
posed for Pro.6oCao.4oMn03 and order in a similar checkerboard configuration associated with the 
CE-type anti-ferromagnetic structure. Therefore, locally the charge difference and structural order- 
ing between the two Mn is appreciable enough to describe the system at 20 K as "partially charge 
ordered" . 



I. INTRODUCTION 

Understanding the relation between the magnetic or- 
dering, electronic ground state, and crystallographic 
structure in mixed valence manganites has preoccupied 
experimentalists and theorists for over 50 years. Today, 
manganites continue to be an interesting topic in solid- 
state physics because of their high degree of functionality, 
strongly correlated electrons, and ability to exhibit fasci- 
nating physical phenomena such as colossal magnetore- 
sistance (CMR)ii In 1955 WoUan and Koeler completed 
their classic neutron powder diffraction (NPD) study of 
the effects of electronic doping on magnetic, structural, 
and transport properties of the manganese perovskite 
Lai_2;Caa;Mn03i^ Since LaMnOa and CaMnOs have Mn 
with formal valences of III (d'' configuration) and IV (d'^ 
configuration) respectively, the system is hole doped as 
La'^^ is substituted by Ca^"*". The half-doped compound 
is an especially interesting composition to study because 
it is on a metal to insulator transition interface in the 
phase diagram. Moreover, one would expect a ratio of 
Mn^+/Mn'^+=1 due to an equal amount of divalent (Ca) 
and trivalent (La) cationsi^ 

A qualitative interpretation of WoUan and Koeler 's 
results proposed by Goodenough predicted that at low 
temperatures the Mn'^+ and Mn'*+ ions in the x — 0.5 
compound would order in alternate (110) planesii The 
structure resulting from Mn^+/Mn'^+ ordering has been 
termed the charge-ordered (CO) structure and its asso- 
ciated magnetic ordering is known as the CE-type anti- 
ferromagnetic (AFM) structure. Moreover, the Mn'^^ oc- 
tahedra are Jahn- Teller distorted and consequently the 
d^ orbitals order in a zigzag pattern; this is described as 
orbital ordering (00) and is coupled with CO in Goodc- 
nough's model. Since the Mn^+ octahedra do not distort, 
their positions are modulated throughout the structure 



with an amplitude that is commensurate with the crystal 
lattice. 

Since the 1990's, Goodenough's interpretation of the 
charge-ordered state has been revised and challenged by 
new experimental studies. The neutron and x-ray pow- 
der diffraction study by Radaelli et alA- suggest that the 
modulation of the Mn'*+ octahedra quasicommensurate 
at low temperature and still confirms Goodenough's pre- 
dictions of charge /orbital ordering (CO/00) 4i However, 
the CO/00 picture has been challenged by x-ray reso- 
nant scattering (RXD), x-ray absorption near-edge spec- 
troscopy (XANES) and single-crystal neutron diffraction 
studies These studies conclude that the local 

structure around the Mn atoms is incompatible with 
Goodenough's Mn-0 bonding model and that the Mn 
atoms do not have integral valence states. In addition, 
one XANES study found that one of the Mn octahedra is 
Jahn- Teller distorted while the other Mn octahedron lost 
their local tetrahedral point group symmetry^ In short, 
the study concludes that both octahedra are distorted 
but in a different manner making the ionic picture of 
CO/00 unrealistic. 

Meanwhile, the study by Daoud-Aladine and 
Rodriguez-Carvajal suggests that the CO/00 
model is inconsistent with single crystal diffrac- 
tion data of Pro.6Cao.4Mn03, a system similar to 
Lai_2;Caa;Mn03i2i2ii& Daoud-Aladine reinterpreted the 
complex magnetic ordering and the modulated structural 
distortion as a consequence of ferromagnetic (FM) Mn- 
Mn dimer ordering. In this model the formal valences 
of the Mn are only fractional (Mn+^-^+VMn+^-S-'^) 
and all the octahedra are elongated similarly so that 
two octahedra connect along their elongated directioni^ 
The FM Mn-Mn dimer in this model is likened to a 
Zener polaron.^ A Zener polaron is defined as a pair of 
Mn'^+ cations that share a hole on a bridging oxygen 
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and in which a Zener double exchange mechanism is 
responsible for the FM coupling of the Mn atoms>ii*i^ 
It is important to note that orbital ordering also occurs 
within the Mn-Mn dimer model. 

Our study aims to use long-range and short-range anal- 
ysis of neutron powder diffraction data to provide fur- 
ther insight into the modulated distortions seen at lower 
temperatures. Rietveld analysis alone will only yield the 
average structure of the material. Here we combine Ri- 
etveld analysis with a local structure analysis technique 
called the pair distribution function (PDF). The PDF is 
the Fourier transform of the total scattering pattern con- 
taining both Bragg as well as diffuse scattering. This ap- 
proach will result in short, medium and long range struc- 
tural information, depending on the refinement range 
Tmax- The PDF approach has been used extensively to 
study the local structure of other manganitesi^'i^'iS'iS'ii 
and a great overview of the technique applied to a variety 
of complex materials can be found here^ 



II. EXPERIMENTAL 

Samples were prepared by conventional ceramic meth- 
ods. Before weighing, La203 (99.99%) was fired at 
1000 °C in O2, Mn02 (99.999 %) was treated in flowing 
oxygen at 425 °C then slow-cooled; CaC03 was dried at 
100 °C. Stoichiometric quantities were mixed and cal- 
cined at 900 °C, then sintered between 1100 °C and 
1380 °C with intermediate grinding and pelletizing. The 
final pellet was annealed in O2 and slow cooled to ensure 
stoichiometry. 

Neutron powder diffraction data were collected on 
NPDF^^ and HIPD neutron diffractometers located at 
the Lujan Neutron Scattering Center at Los Alamos Na- 
tional Laboratory. The sample of about 10 g was sealed 
in a cylindrical vanadium tube with helium exchange gas. 
Data were collected at temperatures between T=300 K 
down to T=20 K in a closed cycle helium refrigerator. 
For PDF analysis the data are corrected for detector 
deadtime and efficiency, background, absorption, multi- 
ple scattering, inelasticity effects and normalized by the 
incident flux and the total sample scattering cross-section 
to yield the total scattering structure fimction, S{Q). 
This is Fourier transformed according to 



TABLE L RT model used in Rietveld analysis of 300 K data. 
Orthorhombic lattice with Pbnm symmetry. Lattice con- 
stants a = 5.4309(1)A, b = 5.4211(1)A, and c = 7.6400(1)A. 
Rwp=4.4% and x^=6.2. 



Atom 


Wyck pos. 


X 


y 


z 


Mn 


4b 


0.5 


0.0 


0.0 


La/Ca 


4c 


-0.0036(2) 


0.0195(1) 


0.25 


01 


4c 


0.0599(1) 


0.4918(2) 


0.25 


02 


8d 


0.7236(1) 


0.2770(1) 


0.0312(1) 



1 " 


300 K 


1 1 
sa& K 


1 

15S K 

■ 71 


1 

100 K 

"A 


1 

so K 




1 

(a) 


i.ee i.9£ 1.86 i.gs i.fls i.es i.ee i.ae i.gs 

d-spBctns [81 

" ! 

.1. t^i .Alt ^. ii ^ . J 


1 1 1 1 1 1 1 1 1 



O.S O.e L.O 1.2 L.4 1.6 l.e 2.0 2.2 2.4 
D-spacing, A 



300 K EES K 155 K 100 K EO K 




^ ^ '^^A-l 3-243-1-)' 3,345.17 3,243.17 3,34 3,1V 3,34 

m d-apacing [A] 



£:-.tji,..j^.:.!:.,-:.:.T.^,,.:!t !^ j.^^,-*- 



0.6 0,8 1.0 1.2 1.4 1.6 1,8 2.0 2,2 2.4 

D- spacing, A 

FIG, 1: (a) NPD pattern at 20 K fitted with the LT-M model. 
Arrows in inset indicate contribution from the room tempera- 
ture phase, (b) The LT-O model fit of the 20 K data. Arrows 
in inset indicate the super-lattice peaks. The bottom line be- 
low the tick marks is the difference between the observed and 
calculated NPD pattern. The temperature evolution in the 
insets are NPD patterns measured from HIPD and the fits 
are done for NPDF data. 



III. RESULTS AND ANALYSIS 



G{r) = - / Q\S{Q) - 1] sin(Qr) dQ. (1) 
The data were terminated at a value of Qmax — 

35 A-i, 

Data processing was carried out using the program 
PDFgetNfSSi Data collection and analysis procedures 
have been described elsewhere«i^ 



A. Rietveld Analysis 

The Rietveld refinements were carried out using the 
GSAS Rietveld codeiSi At room temperature the struc- 
ture of Lao.sCao.sMnOa is single phase and orthorhom- 
bic with space group Pbnm. The parameters of the room 
temperature phase, denoted RT, are listed in TableQ] At 
T=20 K only one phase was observed and found to have 
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TABLE 11: LT-M model used in Rietveld analysis of 20 K 
data with the values of Ax, Ax', and the position of atoms 
refined with P2i/m symmetry constraints. The monoclinic 
cell is based on an orthorhombic cell doubled in the b direc- 
tion. The y and z coordinates for Mn and 02 atoms are fixed. 
Ax = 0.0220(1), Ax' = 0.0103(3), a = 5.47450(7)A, 6 = 
10.88703(13) A, c = 7.51880(8) A, 7 = 89.967(3)°. Rwp=5.7% 
and x^=10.2. 



Atom 


Wyck pos. 


X 


y 


z 


Mn(l) 


2a 


0.0 


0.0 


0.0 


Mr\(9\ 

IV J. 11 \ ^ ) 


2b 


0.0 


0.5 


0.0 


Mn(3) 


4f 


0.5 - Ax 


0.25 


0.0 


La(l) 


2e 


0.4990(5) 


0.0117(4) 


0.25 


La(2) 


2e 


0.4952(5) 


0.5068(4) 


0.25 


La(3) 


2e 


0.0057 - Ax' 


0.2601(4) 


0.25 


La(4) 


2e 


0.0057 + Ax' 


0.7569(3) 


0.25 


01(1) 


2e 


-0.0592(6) 


-0.0133(3) 


0.25 


01(2) 


2e 


-0.0538(5) 


0.4980(3) 


0.25 


01(3) 


2e 


0.5660 - Ax 


0.2532(4) 


0.25 


01(4) 


2e 


0.5660 -1- Ax 


0.7470(4) 


0.25 


02(1) 


4f 


0.2248 - Ax 


0.1440 


0.0337 


02(2) 


4f 


0.7198 - Ax 


0.1151 


-0.0337 


02(3) 


4f 


0.2248 + Ax 


0.6440 


0.0337 


02(4) 


4f 


0.7198 + Ax 


0.6151 


-0.0337 



TABLE IIL Model LT-M' used in Rietveld analysis of 20 K 
data. Atom positions refined with P2i/m symmetry con- 
straints, a = 5.4750(51)A, 6 = 10.8885(1)A, c = 7.5194(1)A, 
7 = 89.950(1)°. Rwp=4.5% and x^=6.5. 



Atom 


Wyck pos. 


X 


y 


z 


Mn(l) 


2a 


0.0 


0.0 


0.0 


Mn(2) 


2b 


0.0 


0.5 


0.0 


Mn(3) 


4f 


0.5157(4) 


0.2469(3) 


0.0078(3) 


La(l) 


2e 


0.0234(5) 


0.2589(3) 


0.25 


La(2) 


2e 


-0.0090(5) 


0.7644(4) 


0.25 


La(3) 


2e 


0.5018(6) 


0.0051(3) 


0.25 


La(4) 


2e 


0.4924(5) 


0.5121(4) 


0.25 


01(1) 


2e 


-0.0577(6) 


0.0003(4) 


0.25 


01(2) 


2e 


-0.0662(5) 


0.4890(3) 


0.25 


01(3) 


2e 


0.5913(4) 


0.2472(2) 


0.25 


01(4) 


2e 


0.5483(5) 


0.7470(4) 


0.25 


02(1) 


4f 


0.2577(4) 


0.1360(3) 


0.0345(4) 


02(2) 


4f 


0.2069(4) 


0.6355(3) 


0.0352(4) 


02(3) 


4f 


0.7372(3) 


0.1129(4) 


-0.0366(4) 


02(4) 


4f 


0.7073(4) 


0.6132(2) 


-0.0294(4) 



antiferromagnetic ordering of the CE-type. In between 
225 K and 100 K there is a mixture of the RT and LT 
phases as shown in Fig. ^ inset. As the temperature is 
lowered, the low temperature phase develops within the 
room temperature phase and its unit cell anisotropically 
distorts. This is similar to other NPD and XRD stud- 
ies noting coexisting phases from around 210 K down to 
low temperatures4*S^ In addition, the room temperature 
phase develops long-range ferromagnetism at Tc=195 K 
and the low temperature phase undergoes an AFM order- 
ing transition at TAr=155 K. Unlike other NPD studies 



TABLE IV: Model LT-O used in Rietveld analysis of 20 K 
data. Atom positions refined with P2\nm symmetry con- 
straints, a = 5.47507(5)A, h = 10.88833(8)A, c = 7.5193(1)A. 
Rwp=4.4% and x^=6.2. 

Atom Wyck pos. x y z 

Mn(l) 4b 0.0176(8) 0.3742(4) 0.2537(6) 

Mn(2) 4b 0.9923(11) 0.8764(4) 0.2445(6) 

La(l) 2a 0.5017(8) 0.3603(3) 0.0 

La(2) 2a 0.5244(8) 0.3838(4) 0.5 

La(3) 2a 0.4923(8) 0.8676(3) 0.0 

La(4) 2a 0.4906(8) 0.8827(3) 0.5 

01(1) 2a 0.0768(10) 0.3692(5) 0.5 

01(2) 2a 0.9448(7) 0.3815(4) 0.0 

01(3) 2a 0.0503(7) 0.8811(3) 0.5 

01(4) 2a 0.9187(9) 0.8852(4) 0.0 

02(1) 4b 0.8020(5) 0.5108(2) 0.2817(4) 

02(2) 4b 0.7265(0) 0.7605(2) 0.2885(4) 

02(3) 4b 0.2289(10) 0.7364(3) 0.2176(5) 

02(4) 4b 0.7482(5) 0.0104(2) 0.2846(3) 



our data does not show any contribution from the rem- 
nant room temperature phase at 20 K (Fig. inset). 

Small satellite peaks observed at 155 K are indicative 
of the orbital ordering transition and cannot be fit with 
model describing the room temperature phase, denoted 
model RT. (Fig. Wp). The most intense satellite peak at 
q = [1, 1/2, 1] is created by the modulated displacement 
of the Mn octahedra. Three structural models were used 
to fit the satellite peaks in the 20 K data. The first model, 
denoted model LT-M, is the structural model describing 
CO/00 at low temperatures as defined by Radaelli et 
al.^ In model LT-M the monoclinic lattice is doubled in 
the b-direction and the bond lengths around the Mn'*+ 
are fixed to be approximately 1.915A. Two modulation 
amplitudes Ax and Ax' are applied to the positions of 
the equatorial O, Mn'*+, and half of the La/Ca atoms 
as shown in Table ITU Along with Ax and Ax', the re- 
mainder of the atoms' positions are refined with P2i/m 
symmetry constraints. The goodness of fit parameters 
obtained were a of 10.2 and an Rwp of 5.7%. 

However, if these constraints are loosened to allow 
atomic positions to adjust freely within the symmetry al- 
lowed by the space group Plxjm but with no additional 
constraints on the Mn'^^-O bonds, the Rwp and val- 
ues improve to values of 4.5% and 6.5 respectively. This 
model, denoted LT-M', is less constrained than model 
LT-M because it does not fix the Mn^^-O bond and does 
not refine modulation vectors Ax and Ax'. Although the 
goodness of fit parameters improve in model LT-M' one 
should note, that the number of parameters refined in the 
non-linear least squares fit arc twice the number as those 
refined in the model LT-M. The parameters of model LT- 
M' are listed in Table ITTll and the Mn-0 bond lengths are 
listed in Table |V| The LT-M' model has many similari- 
ties with the LT-M model except for the distribution of 
Mn^'^-O bond distances. The resulting bond distances 
give some indication that the Mn^+ octahedra may not 
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TABLE V: 


Structural information from Rietveld analysis of models RT, LT-M, LT-M', and LT-O. 


The angles and lengths are 


concerned only with the Mn to equatorial bonds. Only model RT is 


used for 300 K data, the rest for 20 K data. 


Model 


Mn-Oe, 


bond distances (A) 


Mn-Oe, 


bond angles (°) 


RT 


Mn-0 


2 X 1.94562(1) 


Mn-O-Mn 


161.3861(1) 




Mn-O 


2 X 1.94321(1) 






LT-M 


Mn^+(l)-0 


2 X 1.93789(6) 


Mn^+(l)-0-Mn'*+ 


157.130(1) 




Mn='+(l)-0 


2 X 2.09058(6) 




160.1521(9) 




Mn3+(2)-0 


2 X 1.90579(6) 


Mn='+(2)-0-Mn*'+ 


161.2382(8) 




Mn3+(2)-0 


2 X 2.08533(6) 




162.2580(7) 




Mn''+-0 


2 X 1.91551(6) 








Mn''+-0 


2 X 1.91464(6) 






LT-M' 


Mn^+(l)-0 


2 X 1.91209(5) 


Mn^+(l)-0-Mn''+ 


158.7380(7) 




Mn^+(l)-0 


2 X 2.06250(5) 




165.4079(4) 




Mn3+(2)-0 


2 X 1.87984(5) 


Mn=*+(2)-0-Mn''+ 


159.3062(6) 




Mn3+(2)-0 


2 X 2.03272(5) 




162.4260(7) 




Mn*+-0 


1.87014(5) 








Mn^+-0 


2.01323(5) 








Mn'»+-0 


1.92520(5) 








Mn''+-0 


1.95823(5) 






LiL-U 




i.ooooHi ) 


Mn(i j-Kj-NLTiyZ ) 


1 C7 ^i'7AO^^ \ 

iO 1 .0 l4:Z[L) 




Mn(l)-0 


1.99869(1) 


Mn(l)-0-Mn(2) 


162.7741(1) 




Mn(l)-0 


1.91047(1) 


Mn(l)-0-Mn(l) 


161.1720(1) 




Mn(l)-0 


2.01610(1) 


Mn(2)-0-Mn(2) 


163.4689(1) 




Mn(2)-0 


1.95475(1) 








Mn(2)-0 


2.01104(1) 








Mn(2)-0 


1.87888(1) 








Mn(2)-0 


2.00133(1) 







be as isotropic as assumed in the CO/00 state. This 
result anticipates a more rigorous structural analysis in- 
cluding a short-range order technique to give some more 
insight. 

The third model, denoted LT-O, is the structural 
model used to describe the Mn-Mn dimer picture as de- 
fined by Daoud-Aladine et ali^ In model LT-O the or- 
thorhombic lattice is doubled along the 6-axis and the 
atoms positions are refined with P2inm symmetry con- 
straints. Parameters from model LT-O are shown in Ta- 
ble IIVI and resulting structural information in Table 
The third model yields a of 6.2 and Rwp of 4.4%, 
which are approximately equal to those of model LT-M' 
and slightly better than those of model LT-M. Again, the 
number of parameters refined in the least squares fit is ap- 
proximately twice the amount in model LT-O compared 
to the LT-M model. Tab. |V| lists the bond lengths and 
angles of all the models obtained from Rietveld analysis. 

Finally we have a look at the high d-spacing part of 
the diffraction pattern (Fig. |2J|. The magnetic structure 
was fit by an AFM phase of the CE-type using a super 
cell twice the size of the chemical unit cell determined at 
20 K. The symmetry of the magnetic phase is consistent 
with the nuclear phase. As the Shubnikov group used is 
P112[/m which constrains the moments to be in the a-b 
plane. Four different Mn sites were used to create the 
pattern and two different moments were obtained. The 
moment value associated to the so-called Mn'^^ site in the 
CO/00 model is 2.838(11) ^s/Mn and that associated 
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FIG. 2: NPD pattern of 20 K data measured on HIPD low 
angle detector banks. Upper tick marks indicate contribution 
from the CE-type AFM phase and the lower tick marks in- 
dicate the nuclear contribution as defined by model LT-M'. 
The inset shows the CE-type anti-ferromagnetic ordering for 
the Lao.sCao.sMnOa compound at low temperatures. 

to the Mn4+ is 2.599(10) /is/Mn. 



B. Pair Distribution Function Analysis 

The results from the Rietveld analysis make it obvi- 
ous that analyzing the long-range order in the low tem- 
perature phase is insufficient for deciphering the subtle 
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TABLE VI: Structural information from PDF analysis of models RT for 300 K data and LT-M' and LT-0 for 20 K data. The 
angles and lengths are concerned only with the Mn to equatorial O bonds. rniax=15A for all the refinements. 



Model 


Mn-Oe, 


bond distances (A) 


Mn-Oeq bond angles (°) 




RT 


Mn-O 


2 X 1.942(2) 


Mn-O-Mn 


161.35(2) 




Mn-0 


2 X 1.949(2) 






LT-M' 


Mn^(l)-0 


2 X 1.932(4) 


Mn"(l)-0-Mn^ 


159.2(3) 




Mn-^+^lj-O 


2 X 2.051(4) 




166.7(3) 




Mn^+(2)-0 


2 X 1.917(4) 


Mn''+(2)-0-Mn''+ 


157.8(2) 




Mn^+(2)-0 


2 X 2.113(4) 




161.2(2) 




Mn*+-0 


1.948(6) 








Mn*+-0 


1.845(5) 








Mn*+-0 


1.969(6) 








Mn'*+-0 


1.884(6) 






LT-O 


Mn(^ )-C) 

lVJ.il I _L / 


1.89(4) 


IVJ.1.1 1 J. J V_/ IVllll I 


161.1(4) 




Mn(l)-0 


1.901(12) 


Mn(l)-0-Mn(2) 


161.7(5) 




Mn(l)-0 


2.018(9) 


Mn(l)-0-Mn(l) 


158.2(4) 




Mn(l)-0 


1.929(10) 


Mn(2)-0-Mn(2) 


163.5(4) 




Mn(2)-0 


2.08(5) 








Mn(2)-0 


1.88(3) 








Mn(2)-0 


2.01(3) 








Mn(2)-0 


1.95(3) 







changes in the Mn octahedra ordering. Therefore, we an- 
alyzed the information hidden in the diffuse scattering of 
the compound through PDF analysis. The PDF can be 
calculated from a structural model using the relation 



Gcir) 



Anrpo, (2) 



where the sum goes over all pairs of atoms i and j within 
the model crystal separated by r^j . The scattering power 
of atom i is bi and (6) is the average scattering power of 
the sample. To account for the cutoff of S'(Q) at Qmax, the 
calculated function G(r) is then convoluted with a termi- 
nation function, sm(QmaxT/r). Refinements presented in 
this paper were carried out using the program PDFFIT— . 
The program allows to refine structural parameters such 
as lattice parameters, anisotropic atomic displacement 
parameters, position and site occupancies. Even though 
this is similar to the results of a Rietveld refinement, one 
needs to realize that the structural model obtained from 
PDF analysis is strictly only valid for length scales cor- 
responding to the r-range used for the refinement. This 
opens up the possibility to study the local structure on 
different length scales by varying the r range refined. In 
addition to structural parameters, there are two other 
corrections to the calculated PDF: First the finite reso- 
lution, AQ of the instrument leads to a dampening of the 
PDF intensities by exp(— AQ^r^/2). The second correc- 
tion accounts for changes in the PDF peak width. The 
PDF peak width for a pair of atoms i and j is calculated 
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The first term cr^^- is the PDF peak width of the struc- 
tural component due to the atomic displacement param- 
eters. The next two parameters, S and 7, determine the 
sharpening of near neighbor PDF peaks due to corre- 
lated motion, in other words the tendency to move in 
pha,seM^ Finally the parameter a determines the PDF 
peak broadening at very high distances r due to the in- 
strument resolution. Magnetic correlations are neglected 
in the PDF analysis because the magnetic contribution 
fades quickly with increasing Q compared to the nuclear 
contribution. 

Modeling of the PDF data included two objectives: 
First to observe the differences between the 300 K and 
20 K PDF data in order to determine how the local 
structure is affected and second to determine a structural 
model to describe the 20 K data applying a strategy sim- 
ilar to the one used in the Rietveld analysis. The PDF 
data obtained at T=300 K and T=20 K were fit using RT 
model over the range up to rmax^lsA corresponding to 
the distance across 3-4 MnOg octahedra. Note that the 
nearest neighbor peaks in the PDF are negative because 
the neutron scattering length of Mn is negative (see Fig. 

and b). One also observed a large peak at very low 
values of r in the PDF which is purely due to noise prop- 
agating through the Fourier Transform. As expected, 
the data taken at T=300 K are much more accurately 
described by a single-cell, orthorhombic lattice than the 
20 K data. Since the lattice parameters and atom posi- 
tions are refined for both data sets, the differences seen 
in Fig. |3|3 confirm that the satellite peaks found in the 
NPD pattern cannot be modeled by a single cell and give 
rise to local structure distortions affecting the profile of 
the PDF peaks at 20 K. 

After comparing the room temperature and low tem- 
perature data, the low temperature data were refined 
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using two different structural models, all of which used 
the doubled cells obtained from Rietveld analysis as ini- 
tial structures. First, we used model LT-M which is the 
classic CO/00 structural model constructed with space 
group P2i/m, in which only the modulation vectors Ax 
and Ax' were refined. Again the refinement range ex- 
tends to rniax=15A, yielding a high Rwp value of 13.4% 
and values of 0.0208(2) and 0.0165(5) for Aa; and Ax' 
respectively. The second step of the refinement in this 
model loosens the previous modulation vector constraints 
and allows all the atom positions to move with P2i/m 
symmetry constraints. This second step is similar to the 
less constrained LT-M' structural model used in the Ri- 
etveld analysis except that the atom positions were re- 
fined after the modulation vectors Ax and Ax' were ap- 
plied. The second set of refinements yields a better Rwp 
value of 7.3% as shown in Fig. and we refer to this 
model as LT-M' since we allow the atom positions to 
move with P2i/m symmetry constraints. 

The second model is the LT-0 model, which allows the 
atoms to move in any direction but with P2inm symme- 
try constraints. The result is shown in Fig. The 
Rwp=12.8% indicates a much worse agreement of this 
model compared to the LT-M' model when refined up to 
rmax=15A. However, if we extend the refinement range 




5 10 15 



r(A) 

FIG. 3: PDF refinements over the range rj^jn=1.7A and 
rmax=15A: (a) 300 K data using model RT, Rwp=8.8%; (b) 
20 K data using model RT, Rwp=15.2%; (c) 20 K data us- 
ing model LT-M', Rwp=7.3% and (d) 20 K data using model 
LT-O, Rwp=12.8%. 



of the LT-0 model to rniax=45A (approximately 12 Mn 
octahedra long diameter), we then find an agreement of 
R-wp=12.0% for the LT-0 model which is only slightly 
worse compared to Rwp=10.4% obtained for the LT-M' 
model. This is consistent with the Rietveld results yield- 
ing a similar agreement for both models. This means 
that when a PDF refinement is confined to fit only a 
short-range radius of 15A, it preferred the LT-M' model 
to the LT-0 model but for longer correlations, no such 
preference exists. The resulting structures from the rmax 
of IsA fit for both models are shown in Fig. 0] The 
LT-M' model still supports a checkerboard ordering of 
Mn-0 octahedra. Tab. IVII lists the bond lengths and 
angles of all the models obtained from PDF analysis for 

3-11 rniax=15A. 

Based on the refinement of the PDF obtained at 
T=20 K up to rmax^lsA, model LT-M' gives the best de- 
scription of the local structure of Lao.sCao.sMnOa. How- 
ever, the results of the Rietveld refinements gives no such 
clear distinction between the models. The power of the 
PDF technique is the fact that one can probe short as 
well as medium range order by selecting the r-range used 
in the refinement^ In order to investigate the differ- 
ence between PDF and Rietveld results, we refined the 
20 K PDF data out to rmax=45A using the models LT-M' 
(Fig. [SJl) and LT-0 (Fig. \^). The overall agreement is 
R^p=10.4% for LT-M' and Rwp=12.0% for LT-0. How- 
ever, the agreement for 15A< r < 45A indicated in Fig. 
|5l is practically the same. This is consistent with the 
findings from the Rietveld refinements and we basically 
cannot distinguish between LT-M' and LT-0 as model for 
the medium or long range average structure of the sam- 
ple. On the other hand the local structure (r < 15A) is 
quite sensitive to the differences of the two models. From 
the previous refinements with rmax^lSA as well as the 
Rwp-values shown for r < 15A in Fig. [Sjit is obvious that 
locally the better structural description is model LT-M' 
as we have discussed above. The medium range structure 
obtained via PDF or the long range structure obtained 
via Rietveld refinements on the other hand supports both 
models. 




FIG. 4: (a) Resulting structure from the PDF refinement of 
20 K data using model LT-M' for rmax= ISA. (b) Resulting 
structure from the PDF refinement of 20 K data using model 
LT-O for rmax=15A. 
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FIG. 5; PDF refinements over the range r™j„=1.7A and rmax= 45A: (a) 20 K data using model LT-M', Rwp = 10.4% and (b) 
20 K data using model LT-O, Rwp = 12.0%. The Rwp-values indicated on the panels correspond to the Rwp-values for the 
corresponding range in r. 



IV. DISCUSSION 

One of the main criticisms about the CO/00 model 
from X-ray resonance and absorption studies was that the 
valence around the Mn was never found to be fully 4+ or 
3+. Indeed, our study also concludes the same about the 
Mn valence state as calculated from bond valence sums 
following the relation 



exp 



(4) 



where r and h are the so-called bond valence parameters 
and dij is the distance between the central atom i and the 
neighboring atom j. 28,29 ^pj^g ]-,qjj(J valence calculations 
for the two models shown in Fig. 2|are presented in Tab. 

Ivnl 

In the LT-M' model there are 3 Mn sites, where the 
Mn'^+(1) and Mn'^+(2) sites are Jahn- Teller distorted and 



TABLE VII: Bond valences calculated from PDF models LT- 
M' and LT-O refined to rmax=15A. The estimated standard 
deviation on the last digit is given in parenthesis. 



Model 



LT-M' 



LT-O 



Mn site 

Mn^+(1) 

Mn3+(2) 

Mn'*+ 

Mn(l) 

Mn(2) 



Valence 



3.50(4) + 
3.42(4) + 
3.99(4)+ 
3.81(4)+ 
3.58(4)+ 



the Mn^+ is not. It is evident that the bond valence of 
Mn'^+(1) and Mn'^"'"(2) is an appreciable amount away 
from the 3+ and both close to 3.5+ using the bond va- 
lence sums. However, the Mn^"*" bond valence is actually 
very close to 4+. In the LT-O model obtained from PDF 
analysis, there are two different Mn sites, which are all 
Jahn- Teller distorted and therefore do not correspond to 
Mn'^+ and Mn''+ atoms. However, even in this Mn-Mn 
dimer picture, there are two Mn atoms with different 
bond valences as shown in Table IVIII with one close to 
the 3.5+ value and the other close to 4+. One should 
keep in mind, however, that bond valence sum technique 
has sonic shortcomings in describing the bond valence 
states in the mixed-valence manganites. 

Another result similar with the local structure analysis 
from X-ray data but not incompatible with the CO/00 
view is the local symmetry around the Mn'*+ ion in model 
LT-M' obtained from PDF analysis. As demonstrated in 
other studies, our PDF analysis shows that the local envi- 
ronment around the Mn'*+ atom is anisotropic as shown 
in Fig. Instead of four equal Mn-0 bond lengths on the 
equatorial plane of the octahedra, there are two different 
lengths in both directions. Therefore, the point group 
symmetry is lowered from a tetrahedral point symmetry 
of Df^ down to Cs, which only has a mirror plane. This is 
similar to the XANES study, where the point group sym- 
metry of the so-called Mn''+ is not isotropic nor tetra- 
hedral but instead monoclinic— This suggests that the 
Mn^+ is slightly ferroelectric because the Mn^+ is not 
centered within the octahedra. However, no macroscopic 
polarity should be detected because the other Mn''+ in 
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FIG. 6: Local symmetry around two Mn sites, one corre- 
sponding to the Mn^+ and Mn*+ from the LT-M' model ob- 
tained from PDF analysis. The Mn"^"*" octahedron is distorted 
in the expected Jahn- Teller way. The Mn*""" is not isotropic 
as expected in the CO/OO picture but off centered. 

the super-lattice is off-centered in the opposite direction 
and therefore should effectively cancel each other out. 
This is an interesting discovery concerning the Mn'*"'' be- 
cause it demonstrates that it can be anisotropic as shown 
by other short-range structure studies but in a manner 
not contradictory to the CO/00 picture. 

Another argument against the CO/00 picture is the 
similar values of the two Mn moments. In our Rietveld 
analysis the moment values are 2.838(ll)/LtB/Mn and 
2.599(10)^B/Mn for Mn^+ and Mn'*+ respectively. Since 
the moment values are close, this would suggest that eg 
electron is localized over two Mn ions as suggested in 
the Mn-Mn dimer model. However, other similar studies 
showing Mn moments with similar magnitudes conclude 
that the small difference in moment values may have to 
do with incommensurability in the charge-ordered state, 
especially in the half-doped manganites.'*''^*' One of the 
studies concludes that the existence of broadening of cer- 
tain AFM peaks are a result of domain boundaries in 
Mn^"*" sublattice, which causes the moment direction to 
flip across the magnetic domain boundary>4 Therefore, 
the resulting incommensurability of the Mn'^"'" magnetic 
sublattice causes the refined moment value of the Mn'^+ 
ion to average out to a lower value. Indeed, inspection 
of the NPD pattern in Fig. |21 shows that some of the 
AFM peaks are broadened more than others and would 
support the idea of incommensurability in our charge- 
ordered structure. 



V. CONCLUSIONS 

Our neutron powder diffraction study on 
Lao.sCao.sMnOs has been analyzed using Rietveld 
and Pair distribution function analysis in order to 



determine the average and local structure in this 
controversial material. Our analysis shows that the 
ionic CO/00 model and the FM dimer model succeed 
to equally model the low temperature charge ordered 
state from our powder diffraction data using Rietveld 
analysis (long range average structure). Therefore, a 
local structural analysis technique was necessary to find 
the coordination and bond valence state around the 
Mn ion as in X-ray local structure techniques. Our 
PDF analysis goes a step further than analyzing the 
local environment around the Mn atom (up to 4 A). By 
carrying out the PDF refinements to rniax=15A and to 
rmax=45A, we can create a structure that arranges the 
dissimilarly distorted Mn octahedra to bridge the gap 
between local structure analysis and long-range analysis 
from Rietveld. 

The combined Rietveld/PDF analysis shows that the 
local structure of Lao.sCao.sMnOa at 20 K does not sup- 
port the classic, strictly ionic model but does not fit the 
Mn-Mn dimer structure proposed for Pro.6Cao.4Mn03 ei- 
ther. Instead, our PDF data are locally best fit by a 
model with P2i/m symmetry (LT-M') without the con- 
straints necessary to create the ionic, integral valence 
picture. This structural model supports two different 
Mn octahedra that are both distorted but in a different 
manner. One shows an anisotropic distortion associated 
with the Jahn- Teller effect and a bond valence state close 
to 3.5+. The other has monoclinic symmetry, is clearly 
anisotropic, and a bond valence state close to 4-1- . The 
anisotropy in the Mn*+ octahedra implies that the ion 
is slightly shifted in one direction within its octahedron. 
Therefore, while the system is orbital ordered, it is only 
partially charge ordered at low temperatures. These re- 
sults reveal the importance of examining different length 
scales of a complex system and reveal the power of pair 
distribution function when coupled to Rietveld analysis. 
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